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Three-body recombination rates near a Feshbach resonance 
within a two-channel contact interaction model 
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Abstract We calculate the three-body recombination 
rate into a shallow dimer in a gas of cold bosonic atoms 
near a Feshbach resonance using a two-channel contact 
interaction model. The two-channel model naturally des- 
cribes the variation of the scattering length through the 
Feshbach resonance and has a finite effective range. We 
confront the theory with the available experimental data 
and show that the two-channel model is able to quanti- 
tatively describe the existing data. The finite effective 
range leads to a reduction of the scaling factor between 
the recombination minima from the universal value of 
22.7. The reduction is larger for larger effective ranges 
or, correspondingly, for narrower Feshbach resonances. 
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1 Introduction 

Quantum-mechanical three-body systems of identical bo- 
sons exhibit universal features when the two-body scat- 
tering length becomes exceedingly large. In this limit — 
called the universal regime — the properties of the sys- 
tem depend largely on the scattering length alone and 
can be described by a universal one-channel zero-range 
model [10. 

The one-channel zero-range model predicts, in partic- 
ular, that in the limit of large positive scattering length — 
with a shallow dimer — the low-energy recombination rate 
of three identical bosons into a shallow dimer exhibits, 
as function of the scattering length, a geometric scaling: 
characteristic periodic minima in logarithmic scale with 
the period equal to 22.7. 
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Experimentally the recombination rate in the univer- 
sal regime can be investigated in ultra-cold atomic gases 
using the technique of Feshbach resonance which enables 
fine-tuning of the atom-atom scattering length over sev- 
eral orders of magnitude [1, 0; Q • The characteristic scal- 
ing in the recombination rate has been recently observed 
for some alkali atoms, notably '^^K Q and ^Li 0]. 

The one-channel models can only describe Feshbach 
resonances phenomenologically, through a parametriza- 
tion of the scattering length as function of the external 
magnetic field. However recently a two-channel model — 
which provides a natural microscopic description of Fesh- 
bach resonances p, [p — has been suggested for three- 
body calculations [lotll iL Il2j | . In addition, the two-channel 
model has finite effective range — inversely proportional 
to the width of the Feshbach resonance — and can be ex- 
pected to describe the deviations from the universal pre- 
dictions. This model, however, has not been extensively 
compared with experiments. 

The purpose of this investigation is to confront the 
two-channel model with the newly available experimen- 
tal data on the recombination rates, and to determine the 
effects of the finite effective range on the scaling proper- 
ties of the recombination rates. 



2 Two-channel contact interaction model 

Feshbach resonance in scattering of cold atoms arises due 
to an interplay between several coupled channels. The 
essential physics of the resonance can be described by a 
model with two coupled channels [1, @ . 

The system of two identical bosonic atoms in an s- 
wave state is described by a two component wave- function 
Tp. In the following we assume that the energy of the sys- 
tem is below the excitation threshold and refer to the 
ground state and the excited state channels as corre- 
spondingly open and closed. 
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where r is the distance between atoms, and the radial Expanding ([7]) for smaU fcopcn using 
functions Uopen and Wcioscd describe correspondingly the where = 2m*E* /h^, gives the effective range expan- 
open channel where the two atoms are in the ground sion, 
state, and the closed channel where one of the atoms is 
excited. 

In a contact interaction model the radial components 
of the wave-function satisfy two free Schrodinger equa- 
tions, 
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where a and R are the scattering length and effective 
range respectively. 
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where primes denote differentiation with respect to the 
relative distance r, m* is the reduced mass of the two The two-channel contact interaction model thus has a 
atoms, E is the energy of the system, and E* is the negative finite effective range, 
excitation energy from the ground to the excited state 

of the atom. 

The interaction between atoms appear in the contact 
interaction model only through a non-trivial boundary 
condition at small separation between atoms. 
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3 Feshbach resonance and model parameters 

If an external magnetic field is applied to the system, the 
excitation energy E* is substituted with 



E* E* 



(11) 



where the constant /3 parametrises the coupling between where B is the applied magnetic field, and Sfi is the dif- 
thc channels, and Copen and fldoscd are the interaction ference between the magnetic moments of the atom in 



parameters which become scattering lengths in the cor- 
responding channels in the limit of vanishing coupling. 
The model is similar to that of [ll| but is more general 
since we allow Odosed 7^ aopcn- 

For scattering below threshold, < E < E* , the so- 
lutions to the Schrodinger equation ([2]) should be sought 
in the form 



the ground and the excited state. The scattering length 
© is then a function of the magnetic field. 
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where k{B) = y/2m*{E* - S^j.B)/h?. 

The scattering length diverges at a critical value of 
(4a) the magnetic field, Bq, given by 
(4b) 



open COI1S1j3jI11jS , ^opcn 



y^2m*E/h'^ 

and Kciosed = ^/2m*{E* - E)/h?. 

Inserting ansatz (Hj) into the boundary condition ([3]) 
gives 
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The homogeneous linear system ([5]) has non-trivial 
solutions only when the determinant of the matrix van- 
ishes 



Expanding a{B) in the vicinity of Bq gives 
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Solving this equation gives the scattering phase, 
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The characteristic dependence ([15)) of the scattering 
length on the magnetic field is referred to as Feshbach 
resonance Q and is well known empirically. The two- 
channel contact interaction model is thus able to natu- 
rally describe this phenomenon. 
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In the vicinity of the Feshbach resonance the effective 
range is inversely proportional to the width AB of the 
resonance. Indeed from © with a = oo we get 
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in agreement^ with Q . This value holds precisely at the 
resonance only. The i?-dependence of R is, from ^TC 
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Fig. 1 Scattering length a between two Na atoms as a 
function of external mag neti c field strength B with two- 
channel and empirical p5|l fits to experimental data from 
m. The Feshbach resonance is seen at Bo = 907 G with width 
AB = 0.7 G. 



Given AB, Bq and the background scattering length 
Qbg from experiment, the equations (ITil) and (|16p can 
be solved for the model parameters aopon, ^closed and /3, 
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where 



E* ~ S^Bq 
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The value of E* cannot be determined uniquely from 
the equations. It can, however, be found by fitting © to 
experimental data points a{B) as done on figure[TJ How- 
ever, the value of E* does not affect the final observables 

^ Note that m* — m/2 where m is the mass of the atom. 



significantly, provided it is greater than SfiBo and is of 
the order the hyperfine splitting. 

Figure[T]shows experimental data for a Feshbach reso- 
nance in ^■^Na along with fits ([9]) and (jlSp . The phe- 
nomenological expression (fT5|) gives Bq — 907.0 G, AB = 
0.71 G, while the two-channel expression ^ gives Bq = 
907.1 G, AB = 0.69 G when using the expressions (HH) 
and (fT6)) . Choosing E* — 23 /ieV provides a good fit. 
Varying E* has relatively little infiuence on the values 
of Bq and AB provided it lies within this order of mag- 
nitude. 

The phenomenological and the two-channel curves 
are virtually identical, so in the following we use the ex- 
perimental parameters Ubg, AB, Bq and Sfi and deter- 
mine the parameters of the two-channel model through 
the equations (|T^ . 



4 Three-body hyperspherical adiabatic method 

We describe the system of three particles using hyper- 
spherical coordinates defined from the Cartesian coordi- 
nates ri,rj,rk, of particles k as [l3| 
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where {i,j,k} are cyclic permutations of {1,2,3}, p is 
the hyperradius and a is one of the hyperangles, the 
four remaining hyperangles determine the direction of 
Xi and Hi. All five hyperangles are denoted collectively 
as f2. The constant m is a mass-scaling parameter that 
we choose to be the mass of the atoms such that pi — ^ 
for all i. The index i is referred to as the Jacobi-index of 
the chosen coordinate set. 

We expand the wave-function 'F on adiabatic basis 
states ^n(/3, >!7), 

irip, Q) = ^ Q)^ (25) 



where "Pnip, are solutions to the hyperangular equa- 
tion ll 
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where A is the grand angular momentum operator in 
hyperradial coordinates and A„ is the eigenvalue, the in- 
dex n refers to the adiabatic channel. The hyper-radial 
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functions fn{p) satisfy the hyper-radial equations 
A„ + 15/4 
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where the adiabatic couphng terms P and Q are defined 

by 
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with brackets indicating integration over all hyper-angles. 

Contact interaction potentials are equal to zero for 
p > 0, and are specified by applying boundary conditions 
on The boundary condition ^ in hyperspherical co- 
ordinates becomes 
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where Oi 

coupling strength between the open and closed channel of 



lOpeni 



=0 

(29) 

fli, closed are scattering lengths for, and /3i 



particles j and k, and <1> ~ 



•fclosod 



is a two-component 



wave-function that describes the two interaction chan- 
nels. Since the particles are identical, in the following we 
choose a specific Jacobi-set and suppress the indices. 

We use Faddeev decomposition of the hyperangular 
wave- function with s-wave states only, (!> = (j)i + (j32 + 4>i, 
with 

(30) 
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where N and N are amplitudes and where 
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The solutions (jjj and (j)k must be expressed in the same 
coordinate system as in order to use the same bound- 
ary condition. This is done by a rotation and projection 
operator TZ jjj] 
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^>(a,) = (/'(aO-H27^[0fe](a,). (34) 



where ipik = 



Plugging this solution into the boundary condition ([29| 
gives a matrix equation similar to ([S]) of the two-particle 
model. Again a solution exists only when the determi- 
nant equals zero, yielding the eigenvalue equation 
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where k is related to the energy separation E* between 
the channels by E* — fi^K^ /2m*, and the function /; is 
(with I = open, closed) 



sin ( cc— I 
V 6/ 
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(36) 

where x — v when I = open and x = v when I = closed. 
This implicitly defines the eigenvalue v{p) which is the 
central quantity from which all other results follow. 



5 WKB method for recombination rate 



The recombination rate is given by 



(37) 



where n is the particle density and a is the recombina- 
tion coefficient. We calculate the recombination coeffi- 
cient using the WKB method of hidden crossing theory 
[l^ where a is given as 



a = 8(2^)23^3-^ lim:^. 



(38) 



with the wave number k defined by E' = h^k"^ /2m* . The 
transition probability P{k) is 



P{k) =4e~2^sinM 



where 



A + iS= / dpJk^- 
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(40) 



and the integral is taken along a contour in the complex 
p-plane connecting the adiabatic channel n = 1 corre- 
sponding to three free particles to the channel n = 
corresponding to a shallow dimer and a free particle. 
Figure [5] is a visualization of the contour along which 
the integral is calculated. Note that this expression has 
included the Langer correction term 1/4/3^ [l^ which in 
this context modifies the effective radial potential 
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The Q„n terms are not included in hidden crossing the- 
ory 

The integration path must enclose a branch-point ph 
in order to connect channels. The branch-point is found 
by solving [Ij 

for (complex) ui, and evaluating pt — p(yb)- In the one- 
channel model the branch-point is at pb ~ (2.592 + 
2.974i)y/Ia. For the two-channel model p{v) is given 
only implicitly and (j42[) must be solved numerically. The 
above value for pb is, however, approximately correct. 
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Fig. 2 Illustration of the integration path taken in the hid- 
den crossing method. 



Furthermore, in the two-channel model there exist 
a second branch-point between the channel n — and 
n = 1 since the eigenvalue equation (|35p is second order 
in v. This means that the transition probability p9p is 
modified. One has to take the coherent sum of the con- 
tributions from the two paths 



P{k) = \2e-^' sinZ\i 



2e"'^" sinZ\2| 



(43) 



where the indices denote the path taken. However, the 
scattering length has to be rather small compared to 
the eff'ective range for this additional contribution to the 
recombination coefficient to be noticeable. 



6 Thomas effect and recombination rate scaling 



It has been shown [T^ that incorporating a fi- 
nite effective range in the single-channel approach re- 
moves the Thomas effect. Indeed, introducing the effec- 
tive range correction in the boundary condition [ISJ leads 
to the single-channel eigenvalue equation 
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where that imaginary solution — and correspondingly the 
Thomas collapse — is removed at small distances. 

However, in the two-channel contact interaction model 
— despite its finite effective range — the Thomas collapse 
persists. Indeed, at p = the eigenvalue equation ([55]) 
turns into the square of the (corresponding limit of) the 
single-channel equation (H^ , 



7r\ 8 . / 7r\ 
v— ^ sm i^— 



0, 



(48) 



where the imaginary root Vs and hence the Thomas effect 
are still present. Therefore the two-channel model — ^just 
like the unmodified one-channel model — needs a cut-off 
at small distances. 

This imaginary root in the one-channel model leads 
to a characteristic scaling of the recombination rate as 
function of the scattering length a [l|, (TJ] . Indeed, since 
a is the only parameter of the model, the real part of the 
WKB integral (^(11) can, in the large a and low energy 
limit, be estimated as 



In the one-channel model the angular eigenvalue equa- 
tion is given as [H, [O] 



(i^—^ ^ sin fi^-^ = — - sin fi/— V (44) 
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which can be easily obtained from (1551) by setting /3 = 0. 
At small distances, p <C a, this equation has a specific 
imaginary solution 



V, w 1.006i 
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where pc is the regularization cut-off. The recombina- 
tion rate — having sin^ Zi as a factor — will then exhibit, 
as function of a, a characteristic geometric series of min- 
ima at the points determined by the zeros of the sine 
function, 



|i/J In 



n= 1,2,, 



(50) 



This solution creates an exceedingly attractive hyper- The periodic factor of the series is equal 
radial potential at small distances. 



1.262 



(46) 



In the lowest hyper-radial equation (P7)) this leads to 
the collapse of the system, known as the Thomas effect 

[a, in. 



exp 



22.7 



(51) 



In the two-channel model the upper limit in the estimate 
should be modified due to the existence of another 
length parameter, the effective range. Correspondingly, 
the scaling law is also expected to be modified. 
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Table 1 Experimental data for Feshbach resonances for 
three atomic gasses. /is is the Bohr magneton and ao the 
Bohr radius. 



33 



10" 



Experimental Data 
Two-channel model 

,4 




100 



1000 



a/ao 



34 





Bo [G] 


AB [G] 




abg [ao] 


-Rcff [\abg\] 


23Na \£ 


907 


0.70 


3.8 


63 


-21 


i^^Cs 15] 


-11.7 


28.7 


2.3 


1720 


-1.99x10"* 


'''K [6] 


402.4 


-52 


1.5 


-29 


-2.02 


'Li [7] 


736.8 


-192.3 


1.93 


-25 


-3.17 



Fig. 3 Recombination coefBcient a, eq. (|38|l . for ^^Na as a 
function of scattering length a from the two-channel model 
compared with the experimental data from 4^. The theory 
predicts the next minimum to be around ~ lOOOao . 



7 Comparison with experimental data 

We compare the two-channel model with the experimen- 
tal data for cold atomic gasses listed in table [TJ The 
effective ranges are calculated using the formula (TTTl) . 

In figure [3] the result from the two-channel model 
is shown together with the experimental data for ^"^Na. 
The cut-off is fixed by the experimental minimum at 
al = 62ao. The rather large effective range could make 
the finite range effect — reduction of the scaling factor 
down to 15.7 — noticeable. However, at least one addi- 
tional minimum is needed to make a proper comparison. 
Experimental data is not yet available for this range. 

In figure S] the result from the two-channel model is 
shown together with the experimental data for ^^■^Cs. 
The cut-off is fixed by the experimental minimum at 
a* « 210ao. The effective range is very small indeed 
and the results from the two-channel model are virtu- 
ally indistinguishable from the one-channel model with 
the scaling factor of 22.7. The next minimum should be 
found at 03 « 4770ao. 

Figure [5] shows the recombination coefficient for "^^K. 
The recombination minimum = 5650 ± 900 is cho- 
sen to fit the cut-off parameter. The two-channel model 
gives a J = 254ao, with the experimental value of a J = 
(224 ± 7)ao. Overall the two-channel model fits the data 
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Fig. 4 Recombination coefficient a, eq. (|38|l . for ^'^''Cs as a 
function of scattering length a from the two-channel model 
compared with the experimental data from 5]. The theory 
predicts the next minimum to be around 02 ~ 4770ao. 
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Fig. 5 Recombination coefficient a, eq (|38|l . for ^^K as a 
function of scattering length a from the two-channel model 
compared with the experimental data from 6]. The overall 
a* scaling is also shown. 



quite well. Notably the scaling is correct compared to 
experiment. The ratio of minima from the two-channel 
model is 22.2 whereas the experimental value is 25.2±4.1. 
Our result thus lies within the experimental uncertainty. 

The recombination coefficient for ^Li is shown in fig- 
ure IHl The two minima are at = (119 ± ll)ao and 

= (2676 ± 195)ao. The cut-off is fixed by a^, giv- 
ing the two-channel prediction = 125ao. Again the 
theory describes the experimental data very well. The 
two-channel model ratio of minima is 21.4 while the ex- 
perimental value is 22.5 ± 2.6 and again our result lies 
within the experimental uncertainty. 

At small a the experimental recombination coefficient 
is two orders of magnitude higher than the models pre- 
dict (a similar tendency is seen in figure [5]). This is pos- 
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Fig. 6 Recombination coefficient a, eq. (|38p . for ^Li as a 
function of scattering length a from the two-channel model 
compared with the experimental data from [7|. The overall 
a* scaling is also shown. 
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Fig. 7 Recombination coefficient a, eq. (|38p . for '^Li as a 
function of scattering length a from the two-channel model 
with different effective ranges as an illustration of their effect 
on the recombination minima. The scaling factor of ratios of 
minima is reduced the large |-R| becomes. 



sibly due to finite range of the physical potentials which 
cannot be correctly modelled by the presented contact 
interaction potentials. 



8 Finite range effects 

As illustrated on the figures [21 El and [5] the two-channel 
model — with its finite effective range — shows a reduction 
of the scaling factor between the recombination minima 
as compared to the zero-range model value of 22.7 - down 
to 15.7 for 23Na, 22.2 for ^^K, and 21.4 for ^Li. 

A similar reduction from the zero-range result due to 
finite-range effects was also observed in the bound state 
spectrum [21 1. 



For ^^'^Cs there was no noticeable change due to the 
very small effective range. The relatively small reduction 
in the cases of '^^K and ^Li is apparently due to the small 
effective ranges, which in turn is due to the relatively 
large widths of the corresponding Feshbach resonances, 
see table 1. 

In the case of ^^Na the effective range is an order of 
magnitude larger and, correspondingly, the reduction of 
the factor is also larger. Unfortunately at the moment 
there are no experimental data around the second mini- 
mum to confirm the prediction of the two-channel model 
concerning its position of the second minimum. 

On figure [7| we illustrate the effects of the effective 
range on the recombination rate by performing exploratory 
calculations of the rate for ^Li where we vary the effec- 
tive range — by varying AB — keeping the other param- 
eters unchanged. Indeed increasing the effective range 
decreases the scaling factor between the minima. How- 
ever, for a noticeable effect the effective range has to be 
of the order of several dozen ahg which corresponds to 
the resonance width of less than 10 G. 



9 Conclusion 

We have calculated three-body recombination rates near 
Feshbach resonances using the two-channel contact inter- 
action model and compared the results with the available 
experimental data. We have shown that the two-channel 
model is able to quantitatively describe all available data. 

Unlike the ubiquitous one-channel zero-range model 
the two-channel model naturally describes the magnetic 
field dependence of the scattering length through a Fesh- 
bach resonance. In addition it has finite effective range. 
However, despite the finite effective range the Thomas 
effect persists in the two-channel model and a regular- 
ization is still needed. 

The finite effective range — which in this model is 
inversely proportional to the width of the resonance — 
leads to a reduction of the scaling between recombination 
minima as compared to the zero-range scaling factor of 
« 22.7. The effect is, however, relatively small in the in- 
vestigated datasets for '^Li, '^^K and especially for ^'^'^Cs 
because of the relatively small effective ranges or corre- 
spondingly relatively large widths of the corresponding 
Feshbach resonances. 

For ^■^Na data the Feshbach resonance is an order 
of magnitude narrower, and the effective range, corres- 
pondigly, an order of magnitude larger. The two-channel 
model then predicts the scaling factor of 15.4. However, 
this prediction has to be verified experimentally as the 
data in the vicinity of the second minimum is not yet 
available. 
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